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Renormalization group calculations are used to give exact solutions for rigidity percolation on 
hierarchical lattices. Algebraic scaling transformations for a simple example in two dimensions 
produce a transition of second order, with an unstable critical point and associated scaling laws. 
Values are provided for the order parameter exponent j5 — 0.0775 associated with the spanning 
rigid cluster and also for du — 3.533 which is associated with an anomalous lattice dimension 
d and the divergence in the correlation length near the transition. In addition we argue that 
the number of floppy modes F plays the role of a free energy and hence find the exponent a and 
establish hyperscaling. The exact analytical procedures demonstrated on the chosen example readily 
generalize to wider classes of hierarchical lattice. 
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In this letter we re-visit rigidity percolation on a lat- 
tice and show for the first time how renormalization 
group calculations can be exactly performed on partic- 
ular bond-diluted hierarchical lattices in two dimensions 
and show that the transition is second order. This is in 
contrast to the only other exact solution known for the 
rigidity transition, on Cayley tree networks which is first 
order [l|,|l|. 

Phase transitions associated with rigidity have experi- 
mental importance in the elastic behavior in chalcogenide 
glasses [2]) in protein unfolding [^] and in jamming in 
granular materials [5|. Rigidity percolation is similar 
conceptually to the more familiar connectivity percola- 
tion [2i, 2\ J except that instead of demanding a connected 
pathway across the sample, the more stringent condition 
that the connected pathway is also rigid is required. 

Rigidity percolation on networks has been studied 
since 1984 when the concept was first introduced and a 
mean field description proved remarkably accurate [8J, |9J , 
except very close to the phase transition. Subsequent 
work has been mainly numerical [10|, [ll| . The associated 
rigidity phase transition has been most extensively in- 
vestigated on the triangular network in two dimensions 
where numerical studies (using the pebble game algo- 
rithm outlined below) show that the transition is second 
order and described by critical exponents /3 = 0.18 ±0.02 
and di^ — 2.42 ± 0.12 that are distinct from those of con- 
nectivity percolation (/3 = 5/36 — 0.139 and di' — 8/3 — 
2.667). 

Results in three dimensions using the pebble game al- 
gorithm |12| strongly suggest that the rigidity transition 
is first order on a bond diluted face centered cubic lat- 
tice, whereas if angular forces are included whenever two 
adjacent edges are present, the transition is second or- 
der. This is quite different from comrectivity percolation 
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Figure 1. (Color online) An isolated piece of network is used 
shows an isostatic (unstressed) rigid region with edges shown 
by thin lines and an overconstrained (stressed) rigid region 
containing edges shown as thick lines. At the flexible hinges, 
shown as solid circles, angular motion is possible. 



where the transition is always second order in three di- 
mension [6] . Further information comes from Cayley tree 
networks where connectivity percolation is second order, 
whereas rigidity percolation (from a rigid busbar) shows 
a strongly first order transition [2j. 

Some characteristics of networks related to rigidity are 
illustrated in Fig.[T] Such particular network realizations 
are elucidated by exact counting procedures [isl . [l4| such 
as the Maxwell count (used below) and in the pebble 
game algorithm. Both balance constraints against de- 
grees of freedom. The latter finds the rigid clusters and 
the flexible joints between them and also determines re- 
dundant bonds in overconstrained regions, as illustrated 
in Fig. m 

In this letter we show how exact calculatioirs cair be 
performed oir hierarchical networks in two dimensions, 
taking for detailed discussion the Berker lattice (la . Il6| 
shown in Fig. [5] When diluted, this network is one of the 
simplest which captures fundamental generic features for 
rigidity percolation. Each generation is obtained by dec- 
oration of the previous generation, creating an infinite 
sequeirce that can lead to singularities and a phase tran- 
sition. An exact set of equations can be written down, 
relating quantities associated with generations n -I- 1 and 
n, which can be solved at all bond concentrations p by it- 



eration. Most importantly the stable and unstable fixed 
points can be found and the structure of the rigidity 
phase transition can be described by the scaling behavior 
obtained by expanding about the unstable fixed point. 

It is instructive to do a Maxwell count [17] on the first 
three generations of the Berker lattice shown in Fig. [5] 
The number of fioppy modes F is given by the difference 
in the number of degrees of freedom 2V, associated with 
the number of vertices V, and the number of constraints 
which are associated with the number E of edges. How- 
ever, in general not all the edges are independent con- 
straints and so E must be corrected by the number R of 
redundant edges so that 



F ^2V - E + R. 



(1) 



The number of fioppy modes F in Eq. ^ contains the 3 
rigid body motions in two dimensions (two translations 
and one rotation) that become insignificant in the limit of 
a very large number of edges. For the top panel in Fig. [51 
and ignoring R for the moment, i^ — 3 = 2-2 — 1 — 3 = 
while for the second and third panels F — 3 = 2-5 — 8 — 
3 = -1 and F - 3 = 2 - 29 - 64 - 3 = -9 respectively, 
where the negative numbers signify that not all the bonds 
are independent in these rigid diagrams, and we have 
removed the 3 macroscopic fioppy modes on the left hand 
side of the count. Therefore there is a single redundant 
edge in the second panel and 9 redundant edges in the 
third panel of Fig. [2] (one for each of the eight replications 
of the second panel, plus one new one). By removing 
edges randomly from the third panel (i.e. bond diluting), 
first the redundancy is reduced and eventually there is 
no rigid path between the two solid vertices, and rigidity 
is lost. Note that the Maxwell count for Fig. [1] gives 
F — 3 = 2-8 — 13 — 3 = (again removing the three 
fioppy modes) , but as there is one fioppy mode associated 
with the solid vertices, there is also one redundant edge 
associated with the heavier solid edges in the left side of 
the diagram. 

The number of edges E or more precisely bonds iV,^ (p) , 
and the number of vertices V or sites N.^{p) in the n*'* 
generation becomes, for the undiluted case {p — 1) shown 
in Fig. [5] is 



n =0 



n =1 





Figure 2. (Color online) Each generation of the undiluted 
hierarchical lattice is labeled by an index n and n = 0,1,2 
shown here starting at the top. 
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F = iV^(l) = 8", 

y = 7V^(l) = (3.8" + ll)/7. 



(2) 



An important quantity is the mean coordination defined 
by (r) = 2E/V which tends to an asymptotic value (r) = 
14/3 = 4.667 for the undiluted lattice. It is important 
that this quantity be above 4, which is the mean field 
value of the mean coordination needed for rigidity in two 
dimensions [3, \lM- The number of redundant edges is 
(8" — l)/7 so that the fraction of redundant edges for 
large n approaches 1/7= 14.3% in the undiluted lattice. 
For the triangular lattice, this fraction is even higher at 
1/3 = 33.3%. 



Figure 3. (Color online) Showing the four distinct types of 
graphs that lead to a rigid connection between the two solid 
circles. Edges present are shown as solid lines and missing 
edges as dashed lines. Here (A) has all eight bonds present 
and is rigid with one redundant edge and has probability p^, 

(B) has any single edge missing and has probability 8p^{l—p), 

(C) has any pair of edges missing from the three lower (shown) 
or the three upper ones and has probability 6p®(l — p)^ and 

(D) has a triple of edges missing either from the lower or 
upper part of the graph and has probability 2p®(l — p)"^. The 
number of edges in the rigid cluster is indicated by the number 
under each graph. All other graphs (not shown) do not rigidly 
connect the two solid circles. 



For the diluted case an bond is present with probabihty 
p (concentration) and absent with probabihty 1 —p, so the 
probabihty of the two sohd dots being rigidly connected 
in the second panel of Fig.[5]using the weights from Fig. [3] 
is y = / + 8p^(l - p) + 6/(1 - p)2 + 2p^{l - pf = 
2p^ + 2p^ — 3p^. This leads to the relationship between the 
probabilities p„+i,p„ of rigidity percolating in successive 
generations n + 1, n: 



Pn+i = 2pl + 2pl - 3pl 



(3) 



(withpo = p)- The fixed points p* satisfying p„_|_i = Pn = 
p* are the trivial stable fixed points at p* = and p* = 1 
and the non-trivial unstable fixed point p* — 0.9446 — Pc- 
Close to this latter fixed point, Eq. ^ can be linearized 
by differentiating to give {jpn+i —Pc) — •^i(Pn ^Pc) where 
Ai = Vdp\ + 14p^ _ 2^pI = 1.802. 

Using the cluster probabilities and also the number 
of bonds in each rigid spanning cluster from Fig. [31 we 
find from the mean number of bonds that the probabihty 
Pn+i {p) of a bond belong to the percolating rigid cluster 
is given by the recurrence relation 




Figure 4. (Color online) Showing the probability P{p) that 
a bond that is present is also part of the rigid backbone as 
a function of the probability p that an bond is present. The 
four curves shown are for the result of iterating Eq. Q out to 
n — 4,6, 8, 12 terms respectively as the curve steepens. The 
result for n = 12 is shown by the heavier line, and shows 
convergence on the scale of this plot to the singularity at 
Pc = 0.9446. 



Pn- 



-i{p) = -^[^pi 



13p6_i4pJjP„(p) 



(4) 



(with Po{p) = p). Near the unstable fixed point, 
Pn+i{p) = \2Pnip) where A2 = \[5pt + 13pl - Upl] = 
0.9554, showing that the probability of an bond being in 
the percolating cluster renormalizes to zero at the phase 
transition as expected for a second order phase transi- 
tion. From Fig. 3] we can see how the singular behavior 
at the phase transition develops as n increases: n = 12 
appears very close to giving the full singularity. Near 
p = 1, P{p) = p[l - 69(1 -p)^/4 + ...] where the first 
term is just the probability that an bond is present and 
the second term that at least 2 bonds must be removed 
to produce an bond that is present but not part of the 
rigid backbone, as indicated for example in panel (C) of 
Fig. El 

Using the eigenvalues Ai, A2, and A3(= 8) of the lin- 
earized scaling relationships for p, P{p), and iV^ respec- 
tively we obtain exponents v, /3, and fractal dimension- 
ality d, from Ai = b^/"", A2 ^b-^l", and A3 = b'^, where 
h is the dilatation (length scaling) factor between succes- 
sive generations of the hierarchical lattice. However, as is 
typical for such lattices, b is ambiguous |;16,]; so we quote 
only the values of exponents independent of h. These are 
P = 0.0775, which describes how the order parameter 
P{p) goes to zero at the critical point, and the product 
dv = 3.533, which plays a role in hyperscaling which does 
apply here. 

The question of hyperscaling involves the critical ex- 
ponent a that describes the fluctuations associated with 
the specific heat in the system near the phase transition. 
The exponent a is most easily calculated by differenti- 
ating the free energy twice with respect to the the bond 
concentration p, and hyperscaling also relates (when it 



applies) to the free energy. But the question arises as 
to what is an appropriate free energy as rigidity percola- 
tion is not a system described by a Hamiltonian. There 
is strong evidence, outlined below, that the number of 
floppy modes given in Eq. ([T]) serves as the appropriate 
free energy for it. It can be shown that the second deriva- 
tive with respect to p is positive deflnite. For connectivity 
percolation, the free energy can be found as the s — ?> 1 
limit of the s-state Potts model [18] and in that case is 
equivalent to an appropriate version of Eq. ([T]) in which 
redundancy refers to loops or multiple pathways between 
two vertices and the factor 2 is omitted. In this case a 
single floppy mode is associated with an isolated cluster, 
so the free energy is just the total number of isolated 
clusters and of course is an extensive quantity. Finally 
for connectivity and rigidity, these forms of F have been 
used as a free energy for percolation from a busbar onto 
a Cayley tree network [2|. 

Rather than calculate the number F of floppy modes 
directly, it is easier to calculate R in Eq. ([1]) and hence 
determine F . If the number of redundant bonds at gen- 
eration n is Rn+i{p), then 



Rn+l{p)--mn{p)+pl 



(5) 



(with Ro{p) — 0). The factor 8 in Eq. ^ comes from 
the eight fold replication of any redundant bond from the 
previous generation (e.g. going from n = 1 to n = 2 in 
Fig. 12) . The factor p^ comes from additional redundancy 
if all 8 pieces of the graph are rigid (but not necessarily 
redundant). Eq. ^ together with ([2]) provides an iter- 
ative equation for the free energy Fn{p) resulting from 




Figure 5. (Color online) Showing the number of floppy modes 
f{p) (heavy line and left scale ) and also its second deriva- 
tive with respect to p (thin line and right scale) which is the 
specific heat. The vertical dashed line marks the location of 
Pc = 0.9446 which is not at peak of the specific heat. 



F„ip) = 2N:M - N^,{p) + R„{p). 



(6) 



From the eigenvalue Xf{— 8) for the linearized scaling of 
Fn{p) at Pc and large n, we find that the exponent a is 
negative (signifying a cusp) and also establish that the 
hyperscaling relationship to dv is satisfied: 2^ a ~ dv — 
InAi^/lnAi =3.533. 

It is convenient to define the number of floppy modes 
per degree of freedom as fniv) = Fn{p)l['^N^{p)] so that 
< / < 1. Here f{p) is the thermodynamic limit as 
n — > 1 of fn{p)- 111 Fig. [5l we show both / and its second 
derivative with respect to p. Solving Eqs. © at the criti- 
cal point gives f{pc) = 1 — 7pc/Q+pl/Q = 0.00361 and at 
small p, we have f{p) — 1 — 7p/6 + Tp^/AS + ... where the 
term in p^ is the leading correction due to redundancy or 



the onset of dependent constraints. 

In the above treatment of the rigidity percolation prob- 
lem it has only been necessary to consider averages, of 
such things as numbers of stress-carrying bonds, redun- 
dant ones, floppy modes, etc., governed by additive com- 
position rules. Such additivity is absent for processes 
such as percolation conductivity |7| (or elasticity) , where 
probability distributions have to be rescaled. 

For the additive variables of rigidity percolation, prob- 
ability distributions could have been found simply (from 
algebraic recurrence relations for their Laplace trans- 
forms). These can provide further useful information 
e.g. for distinguishing the situations with/without cen- 
tral limit simplicity away from/near the transition. 

To summarize, we have shown for the first time how 
renormalization group procedures can be used to describe 
second order phase transitions involving rigidity percola- 
tion when rigidity percolates on the Berker lattice con- 
sidered here. 

Outstanding questions include more rigorous ap- 
proaches to establish that the number of floppy modes 
F is the appropriate free energy for this problem. In 
addition much insight would be gained by widening the 
scope of the lattices covered. 

In that connection it should be mentioned that the 
Berker lattice discussed here is the simplest member of 
several families for which analytic results have been de- 
rived (all showing continuous transitions) which space 
precludes presenting here. Work continues towards find- 
ing such lattices with a first order rigidity transition, and 
possibly a parameter to tune the transition through a tri- 
critical point. 
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